Light nuclei in supernova envelopes: a quasiparticle gas model 
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We present an equation of state and the composition of low-density supernova matter composed 
of light nuclei with mass number A < 13. We work within the quasiparticle gas model, which 
accounts for bound states with decay time scales larger than the relevant time scale of supernova 
and protoneutron star evolution. The mean-field contribution is included in terms of Skyrme den- 
sity functional. Deuterons, tritons, and '^H(e) nuclei appear in matter in concentrations that are 
substantially higher than those of heavier nuclei. We calculate the critical temperature of deuteron 
condensation in such matter, and demonstrate that the appearance of clusters substantially lowers 
the critical temperature. 



I. INTRODUCTION 

A key ingredient of studies of the formation of neu- 
trino signal in supernova explosions and the supernova 
mechanism itself is the equation of state and composi- 
tion of matter at the densities p < 10^^ g cm~^ and 
within the temperature range < T < 10 MeV [l|. The 
matter at densities below the nuclear saturation density 
is composed of a mixture of nuclei and free (unbound) 
nucleons with a charge neutralizing background of elec- 
trons. The state of the art equations of state that are 
routinely used in the current supernova simulations, in- 
clude the nucleons, the a particles, and a heavy nucleus 
as the independent degrees of freedom 0, [l] . 

Numerical simulations of core collapse supernovas 
demonstrate that about 10^'^ ergs of gravitational bind- 
ing energy is liberated in neutrinos of all flavors [1, |6| . 
Over the time scales relevant for supernovas and pro- 
toneutron stars matter is opaque to neutrinos above den- 
sities lO'^-'^-lO"'^^ g cm""^. The last-scattering surface, 
known as the "neutrinosphere," generates the neutrino 
spectrum of the supernova, which is potentially observ- 
able by the supernova neutrino detectors. The signal 
carries an imprint of physical conditions at the neutri- 
nosphere and can provide information on properties of 
matter under supernova conditions and dense matter in 
general 0, B Neutrino interactions at the neutri- 
nosphere are also of importance for setting the initial 
conditions for possible nucleosynthesis process in super- 
nova winds. 

In this work we focus our attention primarily on the 
composition and the equation of state of dilute isospin 
symmetric and asymmetric nuclear matter. Our goal 
is to introduce a simple setup for treating increasingly 
complex many-body problems related to light nuclei in 
supernova and protoneutron star matter. Our moti- 
vation lies in the Bose-Einstein condensation (BEC) of 
deuterons [13, [IH and a particles 13^ in nuclear mat- 
ter and under supernova conditions |14l |. Furthermore, 
recent computations of two- and three-body binding ener- 
gies in nuclear medium at nonzero temperature and den- 
sity [m allow us a fully quantum mechanical assessment 
of these effects beyond the occupied volume approxima- 



tion. These goals are accomplished by the Quasiparticle 
Gas Model (QGM), described in the following section, 
which treats the nuclei as stable (infinite lifetime) entities 
over the time scales relevant for the supernova and pro- 
toneutron star evolution. Our setup is sufficiently gen- 
eral to allow for resonant states (finite lifetime effects), 
the degeneracy of species, and hence, the possibility of 
boson condensation and medium modifications of bound 
and scattering states in a clustered environment. The 
binding energies of the nuclei ^ < 13 can be computed 
from first principles in free space, e.g., within variational 
theory 15]. In this work we use the experimentally mea- 
sured binding energies of light nuclei [iG*! . Section IIVI 
discusses the modifications of the binding energies of nu- 
clei in the medium on the example of deuterons. Our 
model includes the mandatory mean-field contribution 
to the energy density and the nucleon effective mass due 
to the momentum-dependent mean field (self-energy). 

In recent years a number of studies improved upon 
the equations of state and composition of Refs. [2] and 
0. Instead of using one single heavy nucleus as a rep- 
resentative, an ensemble of nuclei with A < 1000 was 
included in the composition in Refs. [l3l and [l^, treat- 
ing the nuclei as noninteracting Boltzmann gas. These 
statistical ensemble calculations predict nuclei that are 
smaller than those obtained in a single (representative) 
nucleus approximation and they find substantial amounts 
of light nuclei in the composition. The treatment of light 
nuclei has been improved by including the interactions 
among the nucleons and a's on the basis of phase shifts 
(virial expansion) . This approach has been extended 
further to include the three-body bound and scattering 
states ^20(1 and contributions of all nuclei up to A < A [2l[ . 
Neutrino interactions with light nuclei [21| and the ap- 
pearance of A < 4 clusters in a dynamical simulation 
model i22:] have been explored. The changes in the bind- 
ing energies of light nuclei within the Ritz variational 
theory are given in Ref . [2^ . A different view on cluster- 
ization in low-density nuclear matter emerges from the 
studies of the liquid-gas instability, which predicts clus- 
terization of matter into fragments in the spinodal re- 
gion (Ref. [l^l and references therein). The latter process 
is accessible in heavy- ion experiments [25| . 

The implications of the rich and complex composition 
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of matter on the thermodynamics of matter and its ef- 
fect on neutrino transport and other aspects of super- 
nova physics are not yet fully understood. The purpose 
of this work is to advance the study of light clusters in 
supernovas in the following directions. The composition 
of matter is extended to include all stable nuclei up to 
A < 13, a number that is larger than that included in 
studies of light clusters to date. In doing so the quan- 
tum statistics is fully included, i.e., the assumption of 
Boltzmann gas adopted, for example, in nuclear statisti- 
cal ensemble studies is relaxed. It follows then that any 
possible Bose-Einstein condensation (BEC) in clustered 
matter is automatically included in the theory. Indeed, 
we find that there is a BEC of deuterons in the super- 
nova environment. The effect of isospin asymmetry on 
the composition of matter containing nuclei with mass 
number A < 4 is studied. The formalism to address 
this issue is based on expressing the thermodynamical 
potential in terms of a sum over clusters, where each 
term is expressed through the spectral function of the 
corresponding cluster. This method, allows one to ad- 
dress a multitude of effects, such as finite decay width, 
short-lived states, Landau-Pomeranchuk suppression in 
radiation processes, etc. 

To summarize, the novelty of this work lies in the fol- 
lowing: first, the formalism presented here has the advan- 
tage that it represents the contribution of the clusters to 
the thermodynamic potential in terms of their spectral 
functions. Second, we present a complete quantum sta- 
tistical treatment of clusters up to A < 13. Most previ- 
ous work treats clusters as Boltzmann particles. Such an 
approach by default excludes any possible Bose-Einstein 
condensation. Furthermore, the majority of the previous 
works were restricted to clusters up to A < 4, whereas 
we include clusters up to A = 13. Third, we demonstrate 
the Bose-Einstein condensation of deuterons in super- 
nova matter. Fourth, we demonstrate the dependence of 
the A < A cluster abundances on arbitrary isospin asym- 
metry. 

This article is organized as follows. In Sec. HI] we de- 
scribe the quasiparticle model for a mixture of light nuclei 
in symmetric and asymmetric nuclear matter. In Sec. IIIII 
we present the results for the composition, equation of 
state, and deuteron condensation within the QGM. Sec- 
tion |TV] studies the effects of the in-medium modifications 
of the deuteron binding energies and summarizes our re- 
sults. 



II. QUASIPARTICLE GAS MODEL 



of the system into a sum of contributions of clusters 

^{fJ-n,fJ-p,T) ^^na{Ha,T), (1) 
a 

where fin and fip are the chemical potentials of neutrons 
and protons and fia is the chemical potential of a nucleus, 
which is completely characterized by the variable a. The 
chemical equilibrium among the species (baryon number 
and charge conservation) implies that 



fla ^ {A - Z)fXn + Zfip. 



(2) 



At this stage one may either develop a direct pertur- 
bation theory for the thermodynamic potential [26t or 
construct the Green's functions of the theory from ap- 
propriate Martin-Schwinger hierarchy Q and express the 
thermodynamical potential in terms of densities. We fol- 
low the second path. The Martin-Schwinger hierarchy 
is truncated with the help of self-energies such that the 
equation of motion for a nucleus a decouples from oth- 
ers. Then, the thermodynamic potential for each species 
is given by 



naifi,T) = -V dfi'^ na{fi'a,T), 



(3) 



where na{fJ.'^^T) is the number density of a nucleus 
a = {A,Z). By introducing the Fourier transform of 
the finite temperature Green's function iG^{xi,X2) = 
{ilja{xi)ip'^{x2)) , where ip'^{xi) and ipaixi) are the cre- 
ation and annihilation operators of a nucleus a at the 
space-time point xi, we write the densities of species as 



«5a 



dud^p ^^ 
(27r)4 ^ 



dujd'^p 



Sa{(^,p)fF/Bi(^), 



(4) 



where is the degeneracy factor for spin and isospin de- 
grees freedom and the Fermi/Bose distribution functions 
fp/g{uj) account for statistical distribution of species 
with half-integer/integer spin. 



1 ± exp 



OJ 



T 



The spectral function is given by 



Sa{uj,p) = 



[0J-E^{uj,p)Y+n{^.P)l^' 

where the quasiparticle energy is 



(5) 



(6) 



We consider matter composed of unbound nucleons 
and light nuclei with mass numbers A < 13 in thermody- 
namical equilibrium at temperature T and nucleon num- 
ber density n. Each nucleus is characterized by its mass 
number A and charge Z, which we collectively denote by 
a = {A,Z). We expand the thermodynamical potential 



Eaiuj,p) 



2Am 



Ba +ReT,aiuj,p} - fia, (7) 



Ba is the binding energy of the nucleus a = {A,Z), 
is its self-energy, and Ta{uj,p) = 2lmEa{uj,p) 
is the spectral width. We further assume that the nuclei 
under consideration are long-lived on the relevant time 
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scales of supernova evolution, i.e., ra{u!,p) = 0. We 
also assume that the real parts of the self-energies are 
constants independent of momentum and frequency, in 
which case they can be absorbed in the chemical poten- 
tial. Finally, we neglect the effects of medium modifi- 
cation of binding energies; we return to this problem in 
the concluding section. Under these approximations, the 
spectral function is given by 



number density or by the electron fraction Yg = rie/n [see 
Eq. pT|) ]. Although large asymmetries are not realized 
in supernovas, a rapid neutronization process eventually 
equilibrates when the electron fraction reaches ~ 0.05 
in protoneutron stars. The specific conditions prevailing 
in supernova matter, e.g., finite neutrino chemical poten- 
tial, will be considered elsewhere. 



Sa{(^,p) = 27r(5 uj 



P 



2Am 



(8) 



and the energy integral in Eq. ^ is straightforward. 
The defining feature of our model is now transparent - 
the density is the sum of contributions from infinite life- 
time quasiparticles (nuclei) characterized by the value a. 
All relevant thermodynamic quantities can be computed 
from the thermodynamic potential Eq. ([1]); the pressure 
and the entropy are given by 



on 



(9) 



The pressure, entropy and other thermodynamical pa- 
rameters of the electron gas are obtained from the ther- 
modynamic potential 



-geT 



(2^ 



\n[r^{^E,{k)^^i,)\, (10) 



where electron degeneracy factor — 2, the electron 
energy is — y^^+mf , where me is the electron mass, 
and /ie is the chemical potential. The electron density Uf. 
couples to the density of baryonic matter via the charge 
neutrality condition 



A. Density functional 

We start with a brief summary of the Skyrme density 
functional. We assume that the nucleons interact via the 
Skyrme interaction, which is given by 



V{ri,r2) = to5{ri-r2) + 



5{ri - r2)fc^ 



+ k''^5{ri-r2) + t2 k' 5{ri ~ r2) k 
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(12) 



where k = (Vi - V2)/2i and k' = -(Vi - V2)/2i, 
whereby n is the density of nuclear matter. The parame- 
ters tij ^2 and ^3 are determin ed p henomenologically. 
We use the Skill parametrization 27] ; the parameter val- 
ues are to = -1128.75 MeV fm^', ti = 395 MeV fm^ 
t2 = -95 MeV fm^ and ^3 = 1-4 x 10"* MeV fm^. The 
(on-shell) quasiparticle spectrum for nucleons is given by 
ei{p) = /2mi + ReS(ei(p),p) — fjn which we take in the 
quasiparticle approximation, i.e., 



2m* 



(13) 



rie - ^ Zna = 0, 



(11) 



where rie = drie/dfj^e- The thermodynamical potential of 
positrons is obtained upon substituting /ie — —fJ-e- The 
thermodynamical potential of neutrinos of a given flavor 
has the same form as Eq. (jlOp . where the neutrino mass 
and the chemical potential appear instead of the electron 
ones and the neutrino degeneracy factor is = 1. The 
thermodynamical potential of antineutrinos is obtained 
in a similar fashion. 



III. RESULTS 



where z = n,p is the isospin index (n, neutrons; p, pro- 
tons), fi'i ^ fj.^ - Rel^{ei{pF,i),PF,i), where pp^i is the 
Fermi momentum. The effective mass of a nucleon is 
computed from 



n/p 



l + ^(ii+t2) + ^(t2-ii)(l±x) 



(14) 

which we use as a correction for the masses of free and 
bound nucleons. The explicit form of the self-energy is 
immaterial, because its value at the Fermi surface can be 
absorbed in the chemical potential; we drop the prime on 
the chemical potentials hereafter. 



Under supernova conditions the electron fraction in 
matter is fixed and the evolution is nearly adiabatic (con- 
stant entropy). Here, to set the stage, we first explore the 
limit where the matter is isospin symmetric and isother- 
mal. This discussion is followed by a study of a more 
general case of arbitrary isospin asymmetries. Below, the 
isospin asymmetry is characterized either by the asym- 
metry parameter x = {"^n — np)/n, where the 
neutron and proton number densities and n is the total 



B. Isospin symmetric matter (x = 0) 

Our numerical procedure uses the tabulated binding 
energies for nuclei with mass number ^ < 13 with half- 
decay times that are larger than the relevant dynami- 
cal time scales associated with supernova explosions |T6| . 
We first compute the effective masses of nucleons and the 
mean-field from the Skyrme density functional with the 
Skill parametrization. This is followed by a computation 
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FIG. 1: (Color online) Abundances of nuclei Yo, — riajn in 
dilute isospin symmetrical matter composed of nuclei with 
mass numbers A < 13 as a function of matter density, in 
units of nuclear saturation density no — 0.16 fm~'^, at T = 10 
MeV. The abundances of clusters decrease with increasing 
mass number. 
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FIG. 2: (Color online) Dependence of abundances of light 
nuclei < 4 on matter density, in units of no, for two 
temperatures T = 10 MeV (solid lines, blue online) and 6 
MeV (dashed lines, red online). Because of isospin symmetry, 
the abundances of protons and "^He nuclei, which are shown 
by squares and triangles, respectively, are nearly identical to 
those of neutrons and nuclei. 



of the partial densities from Eq. Q with the normaliza- 
tion 71 = X^a ^iid constraint ([2]), which provides us 
the chemical potentials of the species. Finally, we com- 
pute the thermodynamical potential ((!]) from which we 
obtain the pressure and the entropy. The effects of mean- 
field and mass renormalization are small at the relevant 
densities. Furthermore, the results shown below are in- 
sensitive to the choice of Skyrme paramctrization. 

Figure [T] displays the abundances of light nuclei, de- 
fined as 



at constant temperature T = 10 MeV as a function of 
density (in units of nuclear saturation density no = 0.16 
fm~'^). At low densities the matter is dominated by nu- 
cleons with a small (about 10%) admixture of deuterons. 
At intermediate densities the deuteron fraction becomes 
larger than that of the free nucleons; even though the 
population of nuclei becomes more significant, those with 
A> A contribute less than 1% to the total density. Next 
to deuterons '^He and "^H nuclei are the dominant species 
in matter. The a-particle abundance does not exceeds 
0.5% percent at any density. Medium modifications of 
binding energies of nuclei shift the balance between the 
abundances of nucleons and light nuclei in the high den- 
sity part of Fig. [TJ Recall that as n/np — » 1 nuclei dis- 



appear asymptotically, leaving a continuum of nucleons. 
Note that numerically the abundances of neutrons and 
protons in mirror nuclei (obtained by an interchange of 
neutrons and protons) differ slightly because of the dif- 
ferences in their masses and binding energies; however, 
these differences are insignificant on the scales of the fig- 
ure. 

Figured] shows the abundances of dominant species for 
two different temperatures. Reducing the temperature 
from T = 10 to T = 6 MeV increases the abundances of 
light species, such as deuterons, and ^He, while the 
abundance of a particles is suppressed (note that here 
we assume that = const., which is a valid assumption 
only in the low-density limit). 

Figure. [3] shows the chemical potentials of species un- 
der the conditions discussed in Fig.[TJ The relative order- 
ing of the chemical potentials follows from Eq. Be- 
cause ^„ ~ ^p, their absolute value scales as fia ^ -AfJ-n- 
The negative sign of chemical potentials of bosonic (in- 
teger total spin) nuclei implies that these are above their 
critical temperature and density of Bose-Einstein con- 
densation (see Fig. S]). The condition /ia(T) = at 
fixed density is first fulfilled for deuterons (the critical 
temperature of Bose condensation scales as Tc ~ M~^, 
where M is the boson mass). The critical temperature 
of BEC of deuterons as a function of density is shown in 
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FIG. 3: (Color online) Dependence of chemical potentials of 
A < 13 mass number nuclei on matter density, in units of 
no = 0.16 fm"^, at T = 10 MeV. The chemical potentials 
decrease with increasing mass number. 




FIG. 4: (Color online) Dependence of critical temperature of 
Bose-Einstein condensation of deuterons on matter density for 
A < 2 matter (solid line, black online) A < 4 matter (dashed 
line, red online), and yl < 13 matter (squares). 



Fig.[4]in matter consisting of nucleons and deuterons only 
(A <2), nuclei with mass number A < A and nuclei with 
mass number ^4 < 13. It is seen that the presence of nu- 
clei with A > 2 reduces the critical temperature of BEC, 
whereby the effect of adding the nuclei with mass num- 
ber 4 < A < 13 has little effect, because their fraction is 
small. Physically, the presence of heavier nuclei reduces 
the fraction of deuterons in matter, i.e., the "effective" 
density of deuterons. Therefore, the critical temperature 
Tc ~ nl^^ for BEC is reduced. 

Figure [5] displays the pressure and entropy as a func- 
tion of density for several constant temperatures. It is 
seen that the pressure is large for large temperatures and 
increases linearly with density. The entropy is largest at 
low densities and decreases rapidly as the density is in- 
creased. It is seen that it scales linearly with tempera- 
ture, as one would expect for degenerate fermionic mat- 
ter. A comparison of our equation of state with those of 
Refs. [2] and [3,] shows that the differences are small, as 
expected. 



C. Isospin asymmetric matter (x 7^ 0) 

In this subsection we study light nuclei in isospin asym- 
metric nuclear matter. We consider proton-deficient mat- 
ter, i.e., < X = [i^n — np)/n < 1, or in terms of electron 
fraction < Fe < 0.5, which is the relevant case in super- 
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FIG. 5: (Color online) Dependence of pressure (upper panel) 
and entropy (lower panel) on the density of matter for tem- 
peratures T = 10 MeV (solid line, black online), T = 8 MeV 
(dashed line, red online), and T = 6 MeV (dashed-dotted line, 
blue online), fcs is the Boltzmann constant. 



novas and neutron stars. The dependence of the abun- 
dances of light nuclei on the electron fraction Yg at fixed 
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FIG. 6: (Color online) Dependence of abundances of nuclei in 
matter with A < 4 on electron fraction Ye for fixed density 
n = 0.041 fm~^ and two temperatures T = 10 MeV (solid 
lines, blue online) and T = & MeV (dashed lines, red online). 
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Ye 



FIG. 7: (Color online) Dependence of the chemical potentials 
of nuclei in matter with A < 4 on electron fraction Fe for 
fixed density n = 0.041 fm~^ and temperature T = 10 MeV 
(solid lines, blue online). The deuteron chemical potential is 
shown also at T = 6 MeV (dashed line, red online). 



density n = 0.041 fm^^ and two temperatures T = 10 
MeV and T = 6 MeV is shown in Fig. [51 Consider first 
deuterons (the arguments below apply equally to a parti- 
cles and other nuclei with equal numbers of protons and 
neutrons). Their abundance is maximal for Y^. = 0.5. In- 
creasing asymmetry reduces the number of protons that 
are available for building a deuteron; consequently the 
number of deuterons reduces with increasing asymme- 
try and in the limit Ye = they are extinct. Asymmetry 
breaks the degeneracy between the abundances of '^H and 
•^He; the abundance of '^He, which requires two protons 
per neutron, decreases most rapidly. The abundance of 
triton ('^H) is nonmonotonic: it first increases because ex- 
cess neutrons can be easily accommodated in nuclei and 
then decreases because the number of available protons 
vanishes. These two effects make a compromise when 
Ye — 0.25, where triton abundance is maximal. Note that 
the ratio of abundances of deuterons to tritons is inverted 
for large asymmetries; indeed, in symmetric nuclear mat- 
ter the deuterons are the second most abundant species, 
while in asymmetric matter their abundances fall below 
those of tritons for large enough asymmetries. Lower 
temperatures are seen to increase the proton depletion, 
•^He and deuteron abundances decrease faster, and the 
increase in triton abundance at Y,, < 0.5 is more pro- 
nounced. 

The dependence of chemical potentials of light nuclei 
on electron fraction at fixed density n = 0.041 fm~'^ and 




I 1 1 < 1 < 1 < 1 1 1 

0.1 0.2 0.3 0.4 0.5 

Yp 

FIG. 8: (Color online) Dependence of pressure (upper panel) 
and entropy (lower panel) on electron fraction Ve at density 
n = 4.1 X 10^^ fm^"* and several temperatures. 

temperature T = 10 MeV is shown in Fig. [T] The be- 
havior of chemical potentials is understood in analogy 
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FIG. 9: (Color online) Dependence of nucleon (n and p) and 
deuteron ('^H) abundances on density of matter at T = 10 
MeV. Results with free space binding energies are shown 
by solid lines (black online), those with medium modified 
deuteron binding energies by dashed lines (blue online). 

to the behavior of abundances discussed above: the nu- 
clei with equal neutron to proton ratio, as well as those 
that require proton excess are disfavored by asymme- 
try and their chemical potentials are negative and large. 
The chemical potentials of tritons (■^H) are nonmonotonic 
functions of Yg because excess neutrons are responsible 
for its increase for Ye < 0.5, while the proton extinction 
for Fe > is responsible for its decrease. 

The dependence of pressure and the entropy on the 
electron fraction is shown in Fig. [51 It is seen that the 
pressure is lowest in the symmetric case and increases 
with the asymmetry; like in the symmetric case larger 
temperatures sustain larger pressures and entropies. The 
entropy increases with asymmetry starting from the neu- 
tron matter limit Yg = 0, an increase associated with the 
onset of new degrees of freedom (nuclei), which is fol- 
lowed by a decrease as one approaches the isospin sym- 
metric limit. 



IV. SUMMARY AND OUTLOOK 

In this article we set up a quasiparticle gas formalism 
to compute the equation of state and composition of di- 
lute isospin symmetric and asymmetric nuclear matter 
for applications to supernova physics. Our key finding 



is that matter is dominated by the light nuclei, such as 
deuterons, tritons (^H), and '^He isotopes of helium. The 
a-particles (^He) contribute less than 1% to the num- 
ber density. Furthermore, we find that in a large portion 
of the density and temperature diagram deuterons form 
a Bose-Einstein condensate. The effect of heavier clus- 
ters is to reduce the critical temperature of Bose-Einstein 
condensation of deuterons. A novel feature of isospin 
asymmetric matter is the enhancement of the abundances 
of neutron-rich nuclei and the corresponding suppression 
of proton-rich ones. This is clearly manifest in the en- 
hancement of triton abundances with increasing asym- 
metry, which makes tritons the most abundant species 
after neutrons in asymmetric nuclear matter. Compared 
to isospin symmetric matter the relative abundance of 
deuterons and tritons is inverted in strongly asymmetric 
matter. 

The present setup is a useful platform for further exten- 
sions of the theory, which we would like to discuss briefly. 
The binding energies of light nuclei are generally func- 
tions of density and temperature. At high densities and 
low temperatures the binding energies are reduced and 
at some critical values of these parameters bound states 
are dissolved (see Ref. [23| and references therein). Thus, 
for example, nuclei will disappear in matter at high den- 
sities leaving behind a uniform nuclear fluid. The critical 
extinction line for deuterons and tritons in the phase dia- 
gram of symmetric nuclear matter was obtained recently 
in Ref. [lll|. In Fig. [9] we show the effect of incorporat- 
ing the temperature-density-dependent binding energies 
of deuterons, computed in Ref. [llj, on the composition 
of matter with mass numbers A < 2. It is seen that the 
high-density asymptotic state of abundances is inverted; 
the abundances of deuterons are larger than the nucle- 
onic abundances for constant, free space, binding ener- 
gies. However, their relative ratios are inverted when 
the reduction of the deuteron binding energies at large 
densities is taken into account. Matter effects will af- 
fect the abundances of other light nuclei in a similar way, 
which will guarantee that the high-density asymptotic 
state corresponding to the continuum of nuclear ffuid at 
saturation density is recovered. 

Apart from the statistical effect of suppression of 
bound state energies in matter further aspects that 
should be incorporated in the model include (i) leptons 
and electromagnetic forces (screening of nuclear charge), 
(ii) onset of (3 equilibrium during the late time dynamics 
of supernovas, (iii) elastic scattering among the light nu- 
clei themselves and with nucleons, and (iv) reactions. Of 
course, larger numbers of nuclei and resonances (nuclei 
with short-decay times scales) can be easily incorporated 
within our model. 

This work was in part supported by the Deutsche 
Forschungsgemeinschaft (Grant SE 1836/1-1). 
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